Comment on "Ruling out chaos in compact binary systems" 
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In a recent Letter, Schnittman and Rasio jl| argue that they have ruled out chaos in compact binary 
systems since they hnd no positive Lyapunov exponents. In stark contrast, we find that the chaos 
discovered in the original paper under discussion, J.Levin, PRL, 84 3515 (2000), is confirmed by 
the presence of positive Lyapunov exponents. 
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Intuitively, black hole binaries are obvious candidates 
for chaotic dynamics. Firstly, they are highly non-linear 
systems and chaos is an expression of extreme non- 
linearity. Secondly, there are isolated unstable orbits 
around a Schwarzschild black hole, and unstable orbits 
are a red flag for the onset of chaos. Surely the orbits only 
become more unstable and more complex when there are 
two black holes. 

There is no question that there is chaos in black hole 
binaries for some range of parameters [^,^ . Yet the au- 
thors of Ref. claim to rule out chaos by finding no 
positive Lyapunov exponents along the fractal of Ref. . 
Before carrying out a detailed calculation, there is rea- 
son to be suspicious of the absence of positive Lyapunov 
exponents. Even the Schwarzschild solution has posi- 
tive Lyapunov exponents for unstable circular orbits . 
Though completely integrable and so non-chaotic, the or- 
bits are unstable. As long as there are unstable circular 
orbits in the PN equations, there are positive Lyapunov 
exponents for these orbits. This is not to say that the 
orbits will definitely be chaotic, just that there will be 
instability and a positive Lyapunov exponent. 

Chaos was detected in the Post-Newtonian (PN) ex- 
pansion of the two-body problem for rapidly spinning 
bodies |^ using the method of fractal basin boundaries. 
One fully expects positive Lyapunov exponents at least 
very near the fractal basin boundaries, if not elsewhere 
in the phase space. We sampled some orbits near the 
boundary and find positive Lyapunov exponents as illus- 
trated by the positive slope in Fig. |l|. The reason for the 
discrepancy seems to be that the authors of Ref. 0] de- 
fine the maximum exponent as "the Cartesian distance 
between the dimensionless 12-component coordinate vec- 
tors [r, f, Si, S2] and [f, r', S'^, 83] of two nearby trajec- 
tories" However, this is not a Lyapunov exponent. 
The Lyapunov exponent is obtained by linearizing the 
equations of motion about a given trajectory. An ap- 
proximation to the Lyapunov exponent can be made us- 
ing the Cartesian distance between two trajectories by 
continually rescaling the shadow trajectory so that it is 
always infinitesimally close to the original trajectory as 
was done by |H . However, the result does depend on the 
rescaling and can give false answers. A more thorough 
treatment of the Lyapunov exponents and their interpre- 
tation is given elsewhere 13] . 


FIG. 1. Determining the Lyapunov exponent for an orbit 
taken from Fig. 3 of Ref. [3] . 

Chaos has not been ruled out rather it has been con- 
firmed by positive Lyapunov exponents. In fact, an even 
stronger claim can be made. The two-body problem has 
not been solved in relativity and it would not be too dar- 
ing to conjecture that it will never be solved. The system 
shows every indication of being fully non-integrable. 

Regardless of how we'd like nature to behave to make 
our lives as observers easier, she may not be so accom- 
modating. Understandably there is a sense of urgency 
to the removal of any obstacles to the detection of grav- 
itational waves. However, chaos need not be a terrible 
obstacle. There are ways that chaos can enhance signals 
or amplify salient features in the sky |^ to aid our quest 
to observe gravitational waves. 
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